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Microscopic model of (CuCl)LaNb207: 
coupled spin dimers replace a frustrated square lattice 

Alexander A. TsirlirO and Helge RosneiO] 
Max Planck Institute for Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 

We present a microscopic model of the spin-gap quantum magnet (CuCl)LaNb207, previously 
suggested as a realization of the spin-^ frustrated square lattice. Taking advantage of the precise 
atomic positions from recent crystal structure refinement, we evaluate individual exchange integrals 
and construct a minimum model that naturally explains all the available experimental data. Surpris- 
ingly, the deviation from tetragonal symmetry leads to the formation of spin dimers between fourth 
neighbors due to a Cu-Cl-Cl-Cu pathway with an antiferromagnetic exchange J4 ~ 25 K. The 
total interdimer exchange amounts to 12 — 15 K. Our model is in agreement with inelastic neutron 
scattering results and is further confirmed by quantum Monte Carlo simulations of the magnetic 
susceptibility and the high-field magnetization. We establish (CuCl)LaNb207 as a non-frustrated 
system of coupled spin dimers with predominant antiferromagnetic interactions and provide a gen- 
eral perspective for related materials with unusual low-temperature magnetic properties. 
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The spin liquid ground state is one of the great chal- 
lenges in condensed matter physics jii^ While all spin liq- 
uids share the absence of the long-range ordering down 
to zero temperature, the details of their magnetic be- 
havior depend on the specific type of spin correlations, 
largely determined by the lattice topology. A spin liquid 
state is readily achieved in many spin-i gapped mod- 
els (spin dimer, alternating chain, two- leg ladder), where 
the gap in the excitation spectrum results from a sin- 
glet ground state without long-range ordering. Two- 
dimensional (2D) frustrated spin systems show more ex- 
otic spin-liquid regimeSj^ but only a limited range of 
model materials has been studied so fari^i^ 

The (CuCl)LaNb2 07 compound is commonly referred 
as an experimental realization of the spin- 2 frustrated 
square lattice (FSL) model. This model entails com- 
peting nearest- neighbor (Ji) and next- nearest- neighbor 
(J2) couplings on the square lattice and shows the spin 
liquid ground state in a narrow range of parameters 
(J2/J1 ~ 0.5)/- The initially proposed tetragonal sym- 
metry of (CuCl)LaNb2 07 exhibits the square-lattice ar- 
rangement of spin-i Cu"'"^ cationsj^ The spin-gap behav- 
ior gave rise to a common belief that (CuCl)LaNb2 07 
is the first experimental example of a spin liquid regime 
within the FSL model. However, the magnetization data 
fit poorly to theoretical predictions.^'^ Moreover, an in- 
elastic neutron scattering experiment evidenced an un- 
usual position of the gap excitation that could be for- 
mally assigned to a dimer with the length of 8.8 A (com- 
pare to 3.8 A and 5.5 A for Ji and J2, respectively).- 
Nuclear magnetic resonance experiments clearly showed 
a lack of tetragonal symmetry, thus invalidating the FSL- 
based descriptioufS. Until recently, several competing pro- 
posals for the crystal structure and the magnetic model 
were available)^"™ but neither of them could fully ex- 
plain the spin-gap behavior and provide a quantitative 
interpretation of the experimental data. 

Our recent structure refinement, based on a high- 
resolution x-ray diffraction experiment^ yielded accu- 



rate and reliable atomic positions for (CuCl)LaNb2 07. 
Computational approache o^^'^^ seem to converge to the 
same structural model, although the calculated inter- 
atomic distances and angles are slightly different due to 
the inevitable shortcomings of density functional theory 
(DFT), especially in strongly correlated electronic sys- 
tems. In the following, we will use the accurate exper- 
imental structural information to derive individual ex- 
change couplings, to establish the microscopic model, and 
to resolve the long-standing puzzle of (CuCl)LaNb207. 
We also consider the transferability of our model to 
closely related materials^-'J^^ that reveal similar inter- 
pretation problems, although their magnetic behavior is 
strikingly different. 

The evaluation of individual exchange couplings is 
based on scalar-relativistic DFT band structure calcu- 
lations within the local density approximation (LDA)i^ 
and local spin density approximation (LSDA)-l-f/ ap- 
proaches. We used the FPLO code with the basis set of 
atomic- like local orbitals.^^ The on-site Coulomb repul- 
sion parameter U^d was varied in a range of 3.5 — 9.5 eV, 
while the exchange parameter Js^ was fixed at 1 eV. 

The exchange couplings are calculated via two comple- 
mentary procedures: 

i) (model approach) The LDA band structure is 
mapped onto a tight-binding (TB) model and further 
onto a Hubbard model in the strongly correlated regime 
ti <^ C/eff, where ti is a hopping of the TB model and 
C/off is the effective on-site Coulomb repulsion in the 
Cu Zd bands (in general, different from U^d applied to 
the atomic id orbitals). In the half- filling regime, the 
low-lying excitations are described by the Heisenberg 
model with antiferromagnetic (AFM) exchange J^^^ = 
Atf/Ucs- This approach evaluates all the exchange cou- 
plings in the system, yet it does not account for the fer- 
romagnetic (FM) part of the exchange. 

ii) (supercell approach) Total energies for a set of or- 
dered spin configurations from LSDA-)-[/ are mapped 
onto a classical Heisenberg model, thus yielding the total 
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FIG. 1. (Color online) Crystal structure of (CuCl)LaNb207: 
the overall view (left panel) and the [CuCl] "layers" in the 
ab plane (right panel). The couplings are labeled according 
to the Cu-Cu distances: the subscript 1 denotes the inter- 
actions between first neighbors, the subscript 2 denotes the 
interactions between second neighbors, etc. 



exchanges Ji. 

The crystal structure of (CuCl)LaNb2 07 is shown in 
Fig. [TJ The [CuCl] layers in the ab plane were initially 
described within a four-fold symmetry that would lead 
to a square lattice of the Cu+^ cations. However, the 
precise structure determination splits these "layers" into 
chains of corner-sharing CUO2CI2 plaquettesJ^ The LDA 
valence band structure^^ is typical for cuprates. Set- 
ting the Fermi level to zero energy, we find: i) the fully 
filled valence bands below —0.9 eV; ii) the half-filled Cu 
3dx2-y2 bands between —0.3 eV and 0.3 eV (Fig. [2]); 
iii) the wide Nb id bands above 0.5 eV. The apparent 
metallicity is a well-known shortcoming of LDA, when 
appfied to strongly correlated systems. The LSDA+U 
calculations provide a correction to the missing correla- 
tion energy in a mean-field approximation and lead to an 
insulating energy spectrum. ^^ 

Despite the lack of the tetragonal symmetry, the spa- 
tial arrangement of the Cu atoms is close to the square 
lattice. Therefore, we label individual exchange cou- 
plings according to the Cu-Cu distance (see Fig. [T]): 
Ji, J{, and J[' run between first (nearest) neighbors, J2 
and J2 run between second (next-nearest) neighbors, etc. 
The TB fit of the LDA band structure (Fig. ^ identifies 
the relevant AFM interactions (Table HI. Surprisingly, 
the leading AFM interaction is J^^^ ~ 34 K (between 
fourth neighbors) establishing the Cu-Cu dimers that 
control low-energy magnetic properties, see Fig.[3l Other 
relevant AFM interactions include J^^'^^ ~ 23 K as well 
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and Jf ^^ of about 13 K. Further cou- 
plings in the ab plane are below 5 K and can be neglected 
within a minimum model. The total number of inequiva- 
lent exchange couplings up to fourth neighbors in the ab 
plane amounts to 12. The model approach evaluates all 
of them, thus simplifying the supercell calculations. 

The supercell calculations evaluated: i) all the short- 
range couplings {Ji,J[, J", J2, and J2) due to the pos- 
sible FM contributions; ii) the relevant long-range cou- 
plings (J4, J4, and J±). Since all these couplings are 




FIG. 2. (Color online) Tight-binding fit (thick dark lines) of 
the LDA band structure (thin light lines). 



relatively weak, they are sensitive to the choice of the 
Usd parameter in the LSDA-I-C/ calculations. Neverthe- 
less, the qualitative scenario is robust with respect to the 
computational method and can be reproduced for a wide 
range of U^d- In TablelH we list the exchange integrals for 
two representative U^d values that give reasonable agree- 
ment with the experimental energy scale, established by 
the saturation field ^oHc2 = 30 T (about 40 K)^ and the 
Curie- Weiss temperature 6* = 10 K. The different Usd 
values are required due to the different double-counting 
correction (DCC) schemes of LSDA-I-C/: around-mean- 
field (AMF), which is the default option in FPLO, and 
the fully localized limit (FLL) that mimics typical cal- 
culations in the VASP code. A similar offset of 3 — 4 eV 
for Usd depending on the DCC has been previously ob- 
served in other Cu-containing compoundsi ^^i^^ The gen- 
eralized gradient approximation (GGA) for the exchange- 
correlation potential produces nearly the same results. 
AMF and FLL generally favor AFM and FM couplings, 
respectively. This can be seen from the Curie- Weiss 
temperatures, calculated in a mean-field approximation 
{9 — J ^j ZiJi, where Zi is the coordination number for 
Ji): 6* = 26 K for AMF at Usd = 4.5 eV and 6* = -17 K 
for FLL at Usd = 8.5 eV. 

The LSDA-|-t/ calculations enable the establishment 
of the qualitative microscopic scenario. We find sizable 
AFM interactions J4, J4, and J±. The nearest-neighbor 
coupling Ji is ferromagnetic (FM). In the following, we 
will use the experimental data to quantify the micro- 
scopic model. After the spin lattice and the relevant cou- 
plings are established from DFT, efficient numerical tech- 
niques evaluate the properties of the respective Heisen- 
berg Hamiltonian and enable the direct comparison to 
the experiment. Prior to this comparison, we will make 
additional comments on the structural origin of individ- 
ual exchange couplings in (CuCl)LaNb207. 

The FM nature of Ji can be traced back to the twisted 
configuration of corner-sharing CUO2CI2 plaquettes. The 
neighboring plaquettes lie in different planes, thus induc- 
ing the very low J^^^ ~ 13 K. The weak AFM contri- 
bution along with the Hund's coupling on the CI site^S 
lead to the overall FM interaction, despite the Cu-Cl-Cu 
angle of 107.1° notably exceeds 90°, where FM superex- 



TABLE I. Exchange couplings evaluated using the model and 
supercell approaches (see text). Leading couplings are also 
shown in Fig. [T] The model approach is based on the hopping 
parameters i; that are used to calculate AFM contributions 



to the exchange J^ 



Ui/U^a with [/eff = 4 eV.= 



The 



supercell approach evaluates the total exchange integrals Ji 
for two different implementations of the LSDA+C/ method: 
AMF {U-id = 4.5 eV) and ELL iUza = 8.5 eV), see text for 
details. 
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change is expected. The leading AFM couplings between 
the fourth neighbors originate from the effective Cu- 
Cl-Cl-Cu superexchange pathway with two short Cu-Cl 
bonds pointing towards each otheri^^ Thus, the strong 
fourth-neighbor coupling is only possible along the [120] 
and [120] directions, while the couplings along [210] and 
[210] are negligible. The difference between J4 and J4 
is due to the more curved pathway for J[ (the Cu-Cl-Cl 
angle of 154°) compared to J4 (162°). Finally, the sizable 
coupling J± is caused by the low-lying Nb Ad states that 
contribute to the bands near the Fermi level. 

The proposed scenario is highly sensitive to the de- 
tails of the crystal structure. Although relaxation within 
the DFT framework yields the correct crystal symmetry 
and the reasonable structural model, fine features of the 
structure are not properly reproduced. In particular, the 
Cu-Cl-Cu angle for Ji is overestimatedjiSrJ^ This overes- 
timate makes Ji AFM, while J4 is largely overestimated 
compared to the experimental energy scale. Thus, the 
experimental structural information is essential to derive 
the correct spin model of (CuCl)LaNb207. 

The resulting spin lattice of (CuCl)LaNb207 is shown 
in Fig. [21 Its remarkable feature is the lack of the mag- 
netic frustration. The couplings Ji, J4, and J4 tend to es- 
tablish columnar AFM ordering with parallel spins along 
b and antiparallel spins along a. Such an ordering is 
further stabilized by weakly AFM next-nearest-neighbor 
couplings. To test the proposed spin model against the 
experimental data and to quantify the exchange cou- 
plings, we perform quantum Monte-Carlo (QMC) simula- 
tions using the directed loop algorithm in the stochastic 
series expansion representation, as implemented in the 
ALPS simulation package^Si The typical lattice size was 
16 X 16 (1024 sites, four sites per magnetic unit cell) and 
allowed to avoid finite-size effects. Magnetic susceptibil- 
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FIG. 3. (Color online) A sketch of the (CuCl)LaNb207 spin 
lattice in the ab plane. The EM coupling Ji is denoted by the 
dotted line, the AEM couplings J2, J4, and J4 are shown by 
the dashed, thick solid, and thin solid lines, respectively. The 
couplings J± run along the c axis and connect the neighboring 
planes. The shading denotes the spin dimer. The open and 
filled circles show the columnar AFM ordering in the case of 
weak dimerization. 



ity and high-field magnetization data were taken from 
Refs. [13 and [3, respectively. 

Starting from the TB results {J-^^^), we first restrict 
ourselves to the J4 — J[ alternating chain model and fit 
the data with J4 = 25 K and J4/J4 ~ 0.5.^^ However, 
the energy can be transferred between the bonds of the 
lattice, leaving some ambiguity for individual J's. For ex- 
ample, we readily obtained another fit with J4 = 25 K, 
J2/J4 — 0.3, and J4 = (Fig. [H the lower J2 is caused 
by the larger number of the respective bonds). ^'^ One 
can achieve similar fits of the data with an even larger 
number of parameters, but the individual interdimer cou- 
plings remain ambiguous. This implies that the available 
experimental data are insufficient to evaluate fine details 
of the (CuCl)LaNb207 spin lattice. The fits evidence 
the intradimer coupling J4 of about 25 K. The inter- 
dimer coupling amounts to 50 — 60 % of J4 and can be 
distributed among different bonds (J4, J2, Ji, and J±). 
To further characterize the spin lattice, inelastic neutron 
scattering experiments on single crystals are desirable. 
Presently available powder neutron data^ point to an in- 
tradimer distance of 8.8 A in remarkable agreement with 
our model that reveals the dimers on the J4 bond (Cu-Cu 
distance of 8.65 A). 

The proposed spin lattice belongs to the family of cou- 
pled spin dimer models. Since the interdimer couplings 
are non-frustrated, the ground state is determined by the 
ratio of the intradimer and interdimer couplings. If the 
intradimer coupling J4 is sufficiently large, the spin gap is 
opened, as experimentally observed in (CuCl)LaNb207. 

The reduced dimerization will close the spin gap and 
lead to a long-range magnetic ordering. This ordering 
is of the columnar AFM type, because FM Ji along 
with AFM J2,J4, and J4 stabilize the parallel align- 
ment of spins along the b direction (Fig. [3]). The 
columnar AFM ordering has been experimentally ob- 
served in the isostructural (CuBr)LaNb207 (Ref. [ij) 
and (CuCl)LaTa207 (Ref. M)- This resuff demon- 
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FIG. 4. (Color online) Fits of the experimental data with 
J4 = 25 K and J2/J4 = 0.3: temperature dependence of the 
magnetic susceptibility (x) in the primary figure and field 
dependence of the magnetization (M) in the inset. Experi- 
mental data are shown with dots, dark lines are the fits. 



understanding of dimer-based quantum magnets with ex- 
otic magnetic behavior. 

In summary, we have proposed a valid microscopic 
model of (CuCl)LaNb207. We argue that this com- 
pound is a system of spin diniers with non-frustrated in- 
terdimer couplings. The intradimer coupling J4 connects 
fourth neighbors and amounts to 25 K. The interdimer 
couplings comprise 50 — 60 % of J4 and are distributed 
among several bonds of the spin lattice. The model is 
in quantitative agreement with the available experimen- 
tal data and naturally explains the spin-gap behavior of 
(CuCl)LaNb207 as a result of the dimerization. The 
limit of the weak dimerization would lead to the colum- 
nar antiferromagnetic ordering, relevant for isostructural 
compounds. 

We are grateful to Oleg Janson and Artem Abakumov 
for supporting discussions. A.Ts. acknowledges the fi- 
nancial support of Alexander von Humboldt Foundation. 



strates a broader scope of the proposed spin model. It 
can be applied to a range of quantum magnets with 
non-trivial properties. However, the accurate determi- 
nation of model parameters for (CuBr)LaNb207 and 
(CuCl)LaTa207 remains a challenging task and requires 
the precise structure determination along with the in- 
terpretation of the magnetization data. A further chal- 
lenge is the explanation of the -i-magnetization plateau 
in a structurally-related compound (CuBr)Sr2Nb30io-^* 
Such studies are presently underway and will improve our 



Note added: After finalizing our manuscript, an in- 
dependent study of (CuCl)LaNb207 appearedi^ Tassel 
et al., identify (CuCl)LaNb207 as a Shatry-Sutherland 
system with ferromagnetic interdimer couplings. Al- 
though the arrangement of dimers resembles our model, 
we note that the experimental magnetization data could 
not be fitted with exclusively ferromagnetic interdimer 
couplings. There are also no experimental indications of 
the magnetic frustration that is inherent to the Shastry- 
Sutherland model. Further experimental studies are de- 
sirable to resolve the remaining discrepancies. 
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FIG. SI. LDA density of states of (CuCl)LaNb207. The Fermi level is at zero energy. 
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FIG. S2. LSBA+U density of states of (CuCl)LaNb207 calculated with AMF DCC and Uad = 4.5 eV for the lowest-energy 
spin configuration (columnar ordering depicted in Fig. 3 of the paper). The Fermi level is at zero energy. 



